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Abstract: We consider how causality and micro-causality are realised in QED in 
curved spacetime. The photon propagator is found to exhibit novel non-analytic be- 
haviour due to vacuum polarization, which invalidates the Kramers- Kronig dispersion 
relation and calls into question the validity of micro-causality in curved spacetime. 
This non-analyticity is ultimately related to the generic focusing nature of congru- 
ences of geodesies in curved spacetime, as implied by the null energy condition, and 
the existence of conjugate points. These results arise from a calculation of the com- 
plete non-perturbative frequency dependence of the vacuum polarization tensor in 
QED, using novel world-line path integral methods together with the Penrose plane- 
wave limit of spacetime in the neighbourhood of a null geodesic. The refractive 
index of curved spacetime is shown to exhibit superluminal phase velocities, disper- 
sion, absorption (due to 7 e^e~) and bi-refringence, but we demonstrate that 
the wavefront velocity (the high-frequency limit of the phase velocity) is indeed c, 
thereby guaranteeing that causality itself is respected. 



1. Introduction 

The purpose of this letter is to consider how causahty and micro- causahty are reahsed 
in quantum field theory in curved spacetime in the light of the discovery of novel non- 
analytic behaviour in the photon propagator due to vacuum polarization in QED [1]. 

These questions have arisen through the resolution of a long-standing puzzle in 
"quantum gravitational optics" [2,3], viz. how to reconcile the fact that the low- 
frequency phase velocity fph(O) for photons propagating in curved spacetime may 
be superluminal [4] with the requirement of causality^ that the wavefront velocity 
Vwf = Vph{oo) should not exceed c. 

The reason why this is a problem is related to the Kramers- Kronig dispersion 
relation [8,9]. In terms of the refractive index n{uj), where t'ph(co') = 1/Re n{uj) 
(setting c = 1), this is 

Ren(oo) - Ren(O) = / — Imn(cu) . (1.1) 

Provided Im n{uj) > 0, as required by unitarity in the form of the optical theorem, 
this implies Ren(oo) < Re?2(0) and hence Vph{oo) > fph(O). The fundamental 
assumption in the derivation of eq. ( |1 . 1| ) is that n{uj) is analytic in the upper-half 
complex UJ plane, which is generally presented (see, e.g. ref. [10]) in flat spacetime as 
a direct consequence of micro-causality. 

The resolution of this apparent paradox is that the Kramers-Kronig dispersion 
relation does not hold in the form ( |1 . 1| ) . We flnd that n{uj) develops branch-point 
singularities on the imaginary axis in the upper-half plane. The consequent modifl- 
cation of eq. (|1 . 1| ) then allows Ren(oo) > Ren(O) and we flnd by explicit calculation 
of the full non-perturbative frequency dependence of the refractive index that n{oo) 
is indeed equal to 1 and the wavefront velocity itself is f„f = c. Remarkably, we flnd 
that these unusual analyticity properties can be traced very directly to the focusing 
property of null congruences and the existence of conjugate points, which are, in 
turn, a consequence of the null energy condition. Conjugate points are points in 
spacetime that are joined by a family of geodesies: at least in an inflnitesimal sense, 
see flg. (|l]). Such inflnitesimal deformations are associated to zero modes that lead 
directly to the branch points on the imaginary axis. 

^In fact, the question of whether causahty could be maintained in curved spacetime even if the 
wavefront velocity exceeds c is more subtle and involves the general relativistic notion of "stable 
causality" [5,6]. See ref. [7] for a careful discussion. 
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Figure 1: Points p and g on a geodesic 7 are conjugate if they can be joined — at least in an 
infinitesimal sense — by a family of geodesies close to 7. The existence of conjugate points means 
that the classical saddle-point solution with the photon connected at p and g to a degenerate e^e~ 
loop squashed onto the geodesic 7 is deformable into a non-degenerate loop as shown. 

The refractive index in QED is determined by the vacuum polarization tensor 
I{^jj{k). On-shell, and in a basis diagonal with respect to the polarizations ii {i = 

1, 2), the relation is 

riiiu) = 1 + ^U,i{u) . (1.2) 

We have, for the first time, evaluated the complete non-perturbative frequency depen- 
dence of the vacuum polarization Ilij{uj) in QED in curved spacetime by combining 
two powerful techniques: (i) the world-line sigma model, which enables the non- 
perturbative frequency dependence of Ilij{uj) to be calculated using a saddle-point 
expansion about a geometrically motivated classical solution and (ii) the Penrose 
plane- wave limit, which encodes the relevant tidal effects of the spacetime curva- 
ture in the neighbourhood of the original null geodesic traced by the photon in the 
classical theory. 

2. The World-Line Sigma Model, Penrose Limit and Null 
Congruences 



In the world-line formalism for scalar QED, the one-loop vacuum polarization tensor 
is given by^ 

roc jrp pT 

nl/°°' = £/ dr,dr,Z{V:jx{r,)]V^,^\x{r,)]) . (2.1) 

Here, V^,£ja;(r)] are vertex operators for the photon and the expectation value is 
calculated in the 1-dim world-line sigma model involving periodic fields x^{t) = 
x^{t -f T) with an action 

S = [ (^^7^.(x)x'^x'^ - m^) , (2.2) 



^For relevant references on the world-line technique [11, 12] (for a review, see [13]) in curved 
spacetime, see refs. [14-18]. We consider scalar QED for simplicity. The extension to spinor QED 
involves the addition of a further Grassmann variable in the action (2.2), but introduces no new 
conceptual issues. 
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where 5^1/ is the metric of the background spacetime.'^ The expectation value is 
represented to 0{a) by an e"'"e~ loop with insertions of the photon vertex operators 
at Ti and T2, as illustrated in fig. (0). The factor Z is the partition function of the 
world-line sigma model relative to flat space. 




II' 

Figure 2: The e+e" loop x^ {t) with insertions of photon vertex operators at n and T2. 

The form of the vertex operators is determined at this order by the classical 
equations of motion for the gauge field in the geometric optics, or WKB, 

limit. ^ That is, 

= xM^(x) (2.3) 

where 

A^{x) = A{x)i,,{x)e'^''^^^ + ■ ■ ■ , (2.4) 

where A is the scalar amplitude, ii are the polarizations and the phase satisfies 
the eikonal equation 

g^'d.Qd^Q = . (2.5) 

The solution determines a congruence of null geodesies where i^^ = d^Q is the tan- 
gent vector to the null geodesic in the congruence passing through the point x. In 
the particle interpretation, k'^ = ujC.^ is momentum of a photon travelling along the 

^Since there is in general no translational invariance in curved spacetime, the effective action 
depends on a point in spacetime and we choose our coordinates such that this is — 0. In the sigma 
model we deal with the corresponding zero mode in the "string inspired" way by imposing [13-18] 

/ dTa;^(T)=0. 
Jo 

This yields a translationally invariant formalism on the world-line that allows us to fix ri — 0. We 
will then take T2 = £,T, < ^ < 1, and replace the two integrals over n and T2 by a single integral 
over the variable ^. 

^There are 3 dimensionful quantities in this problem: the electron mass m, the frequency uj of 
the photon, and the curvature scale R (of mass dimension 2). We work in the WKB limit uj ^ ^/R 
in a weakly-curved background R <C m^. This leaves the dimensionless ratio uj^R/m^ to define the 
high and low frequency regimes. 
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geodesic. The polarization vectors are orthogonal to ^^ and are parallel trans- 
ported: i ■ D s'^ = 0, while the scalar amplitude satisfies 

log^= = (2.6) 

where the expansion 9 is one of the optical scalars appearing in the Raychoudhuri 
equations. 

In order to evaluate the world-line path integral over x^{t), we will need to 
consider fluctuations about the classical geodesic 7. The first step is to set up 
Fermi normal coordinates (m, v, y*) [19] which are adapted to the null geodesic in 
the sense that 7 is the curve (-u, 0, 0, 0) where u is the affine parameter, v is another 
null coordinate and parameterize the transverse subspace. Now, as explained 
in Section 3, the relevant curvature degrees of freedom needed to describe these 
fiuctuations at leading order in an expansion in Rjvr? are captured in the Penrose 
limit of the spacetime around 7 [20,21]. This follows from an overall Weyl re-scaling 
ds^ — > }^ds^ of the metric obtained by an asymmetric re-scahng of the coordinates, 
{u,v,y^) — > {u, X'^v, Xy^), chosen so that the affine parameter u is invariant. This 
leads, for an arbitrary background spacetime, to the plane wave metric 

ds'^ = 2du dv + hij{u)y' y^ du^ - dy'^ , (2.7) 

where hij{u) is related to the curvature of the original metric g^y by hij = —Ruiuj- 
This is the Penrose limit of the background spacetime in the neighbourhood of 7 in 
Brinkmann coordinates. 




Figure 3: The behaviour of the congruence of null geodesies for the (a) Type I and (b) Type II 
plane wave metrics. 

It is now convenient to divide these spacetimes into two classes, depending on 
the behaviour of the null congruence. Introducing a second optical scalar, the shear 
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0" = y ^Di^^i^'fD'^i'' — 6'^, we can write the Raychoudhuri equations in the form: 

- a) = -(^-a)2- $00 + 1^01 . 
Here, we have introduced the Newman-Penrose notation for the components of the 
Ricci and Weyl tensors: $oo = ^Ruu and |\l/o| = ^{Cuiui — Cu2u'i)- The effect of 
expansion and shear is visuahzed by the effect on a circular cross-section of the 
null congruence as the affine parameter u is varied: the expansion Q gives a uni- 
form expansion whereas the shear a produces a squashing with expansion along one 
transverse axis and compression along the other. The combinations ^ ± o" therefore 
describe the focusing or defocusing of the null rays in the two orthogonal transverse 
axes. Provided the signs of $oo i |^o| remain fixed (as in the symmetric plane wave 
example considered in detail later) we can therefore divide the plane wave metrics 
into two classes (see fig. (|^)). A Type I spacetime, where $oo i l^ol are both posi- 
tive, has focusing in both directions, whereas Type II, where <l>oo=t l^ol have opposite 
signs, has one focusing and one defocusing direction. Note, however, that there is 
no spacetime with both directions simultaneously defocusing, since the null-energy 
condition requires $oo > 0. 

It is clear that provided the geodesies are complete, those in a focusing direction 
will eventually cross. In fact, the existence of conjugate points, as described earlier, 
is generic in spacetimes satisfying the null energy condition [5,22]. The existence of 
conjugate points plays a crucial role in the world-line path integral formalism since, 
as explained below, they imply the existence of zero modes in the partition function 
which ultimately are responsible for the Kramers- Kronig violating singularities in 
the vacuum polarization tensor. 

First, we need the explicit solutions for the geodesic equations in the plane wave 
metric (|2.7| ). With u itself as the affine parameter, these are: 

dlv + 2h,,y'duy^ + \dji,,y'y^ = , dlf + h,,y^ = , (2.9) 

the latter being the Jacobi equation for the geodesic deviation vector y*. The solution 
for V determines the eikonal phase [1] 

e{x) =v- ]flijyY , (2.10) 

where = duEiaE°'j is most simply expressed in terms of a zweibein^ E°'i{u) in 

^The zweibein E°'i{u) and its inverse E^a{u) relate the Brinkniann coordinates used here to 
Rosen coordinates , which are particularly well-suited to describing the null congruence but not 
so simple for evaluating the world- line path integral [1]. A particular geodesic in the congruence 
has = E^aX"' and v — Q ~\- j^ijU^V'' for constant In our conventions i and j in this 2d 

Euclidean subspace are raised and lowered with 5ij. 
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terms of which the curvature is hij = —d^EiaE°-j. The polarization 1-forms and 
scalar amphtude are 



These give all the ingredients necessary for the vertex operators, so we determine 



K,.. [x''{t)] = {i/ - n^^y^u) \l exp lu [v - \n^,yY] • (2.12) 



3. World-line Calculation of the Vacuum Polarization 



We have now assembled all the elements needed to calculate the world-line path 
integral ( p.l|) for the vacuum polarization. The fundamental idea is to evaluate this 
by considering the Gaussian fluctuations about a saddle point given by the classical 
solution of the equations of motion for the action 

S = -T+ — J^ drg^,{x)x'^x^-u;Q[x{0]+u^<d[xm, (3-1) 



including the phase Q{x) of the vertex operators which act as sources. In ( |3.1| ), 
we have re-scaled T — > T/m^, which makes it clear that 1/m^ plays the role of 
a conventional coupling constant. In fact, the effective dimensionless coupling is 
actually the dimensionless ratio R/m^. (We have also re-scaled r Tr so that the 
r integral now runs from to 1.) 

The sources act as impulses which insert world-line momentum at the special 
points x(0) and x(^). In between these points, the classical solution is simply a null 
geodesic path and it is straightforward to see that this is given by f = = while 
u satisfies 

fi = -^«(^-{) + ^«W. (3.2) 



The solution is [1] 



~, , r2^T(l-Or/m2 0<t<^, 

u = u(t) = —uo + < 0.3) 
\2ujTai-r)/m^ ^<t<1, 

where uq = ujTC,{1 — S,)/Tn?- This describes an e^e~ loop squashed down to lie along 
the original photon null geodesic between u = ±Uo as illustrated in fig. (p. 

An intriguing aspect of this is that the affine parameter length of the loop actually 
increases as the frequency uj is increased. Technically, this is simply because for 
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higher frequencies the sources impart larger impulses. However, it means that in the 
high frequency limit, the photon vacuum polarization probes the entire length of the 
null geodesic and becomes sensitive to global aspects of the null congruence. This 
is rather counter-intuitive, since we would naively expect high frequencies to probe 
only local regions of spacetime, and appears to be yet another example of the sort 
of UV-IR mixing phenomenon seen in other contexts involving quantum gravity or 
string theory. 



Another crucial point is that while ( |3.3|) is the only general solution, for specific 
values of T (or equivalently u) there are further solutions. These are due to the 
existence of conjugate points on the null geodesic. For values of T for which v{±uo) = 
y\±uo) = 0, there exists more than one null geodesic path between the points 
r = 0, ^ where the impulses are applied. This gives rise to a continuous set of classical 
solutions, which results in zero modes in the path integral for these specific values 
of T: see fig. (|I|). In turn, these produce the singularities in the partition function 
which are responsible for the violation of the Kramers-Kronig relation. Notice, that 
it is not necessary that these more general geodesies lift from the Penrose limit to 
the full metric. 

The perturbative expansion in R/m? about this solution can be made manifest by 
an appropriate re-scaling of the "fields" x^(r). However, this re-scaling must be done 
in such a way as to leave the classical solution u{t) invariant. But this is precisely 
(see ref. [1] for details) the Penrose re-scaling described above, where A is identified 
with the effective coupling Rjw? . This is why the physics of vacuum polarization is 
captured perturbatively in 0{R/rn?) by the Penrose expansion in 0{\) [1,19] of the 
background spacetime around the photon's null geodesic 7. 



Expanding around the classical solution, in the plane wave metric ( p.7|) , we 
find the Gaussian fluctuations in the transverse directions (the u, v fluctuations are 
identical to flat space) are governed by the action 



•1 



= dr 



'0 



^(y^2 + ^2/i.,(%y) + |^%(%V(5(r-0-5(r))] . (3.4) 



At this point, in order to illustrate the general features of our analysis with a 
simple example, we restrict the plane wave background to the special case where the 
curvature tensor is covariantly constant. This defines a locally symmetric spacetime, 
and the corresponding "symmetric plane wave" has metric ( p^.Tp with 

MV = ^i(z/T + ^2(i/')'- (3.5) 

The curvatures are Ruu = c^i +^^1 C'mImi = — C*m2m2 = |(o"2~^i)- The signs of the 
coefficients determine which class the spacetime belongs to. For Type I (focus/focus). 
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(Ti and (72 are both real, while Type II (focus/defocus) has cxi real and 02 imaginary. 
The case cri,(j2 both imaginary is forbidden by the null energy condition, which 
requires R^u > 0. With this background, the geodesic equations determine the 
zweibein E\ = 5ai cos((TiM + Oj) and VLij = — (5jj(Tj tan((TjM + Oj) (the are arbitrary 
integration constants) . 

The core of the calculation reported in ref. [1] is the evaluation of the partition 
function Z and the two-point function {ii ■ x ij ■ x) for the world-line path integral 
in the symmetric plane wave background. A lengthy calculation described in 
detail in [1] eventually gives the following result for the vacuum polarization: 



l-loop 



-5. 



am 



27r 



00+it 



sinh Pi cosh /3, 



n 

1=1 



sinh f3i cosh (3i 



(3.6) 

where f3i = ujT^{1 — ^)(Ji/m?. In deriving ( |3.6| ), we have performed a Wick rota- 
tion T — i> —iT to leave a convergent integral. The ie prescription deals with the 
singularities on the real axis which arise in the Type II case. 




Figure 4: The 71 = 1 (red) and 71 = 2 (blue) zero modes for ^ = i. The points u = ±uo are 
conjugate points for 7. 

This expansion demonstrates clearly the unconventional analyticity properties 
of the vacuum polarization in curved spacetime. When cxj is real, there are branch 
point singularities on the imaginary T axis at the specific values 



mm n 



^ .... , n = l,2,.... (3.7) 

These singularities arise from zeros of the fluctuation determinant Z and have a nat- 
ural interpretation in terms of zero modes, i.e. non-trivial solutions of the fluctuation 
equations that follow by varying ( |3.4| ) with respect to y\ For the special case ^ = |, 
these zero modes are particularly simple to write down: u = u{t) as in ( |3.3| ) while 

y%T) = sin(2n7rr) . (3.8) 

The n = 1 and n = 2 zero modes are illustrated in fig. (1). (For generic ^, the zero 
modes are more complicated: see ref. [1].) This confirms that the zero modes are 
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associated to geodesies that intersect 7 at both u = ±mo, i-e- the conjugate points 
on the geodesic 7. In terms of the vacuum polarization Ill~^°°^{u!) or the refractive 
index nj(a;), these singularities appear on the imaginary axis in the upper-half plane 
in LJ. They have profound consequences for the Kramers-Kronig dispersion relation. 

In the Type II case, where one of the cTj is imaginary, there are also singularities 
on the real T axis. These singularities can also be understood in terms of zero modes, 
but now with imaginary affine parameter u —>■ iu, and can be thought of as world- 
line instantons. These give rise to an imaginary part for n|j'°°^(co') and the refractive 
index n{uj), which corresponds to the tunnelling process 7 —>■ e^e" in the background 
gravitational field. 



4. Refractive Index in Curved Spacetime 



The refractive index follows immediately from the result ( |3.(j| ) for the vacuum polar- 
ization. We present the results for the Type I and Type II symmetric plane wave 
backgrounds separately. 

Type I: This includes the special case of conformally flat spacetimes, where cxi = 
(72 = \/R. Evaluating (|3.6|) numerically gives the result shown in fig. (|^). Here, 
both polarizations are superluminal at low frequencies, with the refractive index 
rising monotonically to the high-frequency limit n{uj) —>■ 1. The wavefront velocity, 
Vwf = Vph{oo) is therefore c, in accordance with our expectations from causality. The 
integrand in ( ^.61 ) is regular on the real axis and so Im n{uj) is vanishing and there 
is no pair creation. 

We can find explicit analytic expressions for n{uj) in these limits [1]. For low 
frequencies, 

aR vl 71 iu'^R 428 /uj^R\2 1 ^ 

ni(uj) = l — r • (4.1) 

^ ^ 27rm2Ll8 14175 189189 V / J ^ ^ 

This series is divergent but alternating and this is correlated with the fact that it 
is Borel summable, with the sum being defined by the convergent integral in (|3.6| ) 
which has no singularities on the real axis. In the high-frequency limit, we find that 
the refractive index approaches 1 from below, with a 1/a; dependence: 



where = (i + ^) for both i = l,2. 
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n-1 Re n-1 




Figure 5: The behaviour oiKeni{uj) — 1 in units of Q:i?/(27rm^), as a function of ^logoj^R/m'^ 
for (a) Type I conformaUy flat case {ni ~ and (b) Type II Ricci flat case [i = 1 red, i — 2 
green). Note that in both cases ni{u}) approaches 1 from below as a; oo. 



Im n_l Im n_2 




Figure 6: The behaviour of (a) lmni{uj) and (b) lmn2{oj) in units of ai?/(27rm^), as a function 
of ^ logijj^i?/m^ for the Type II Ricci flat case. 



Type II: In this case, the integrand in ( |3.6|) has branch point singularities on both 
the real and imaginary axes. The refractive indices therefore have both a real and 
imaginary part. These are shown figs. (§) and (§) for the special case of a Ricci 
flat background with cxi = ia2 = \fR. Notice that the subluminal polarization state 
n2{uj) has the same general form as the refractive index for a conventional absorptive 
optical medium [1,3]. 

In this case the propagation displays gravitational bi-refringence [4, 23] since the 
two polarizations have different phase velocities. In general, the low-frequency limit 
of the refractive index is 

a 1 

n,{uj) = 1 - — ^(lO<l>oo T 4|^o|) + Oicu^) . (4.3) 

for 2 = 1, 2. For Type II, at low frequencies, we find has 

ni,2[uJ) - T 2^^2 T 28350 ^ 85995 ^"'\- 

for the superluminal and subluminal polarizations respectively. In both cases the 
Borel transforms have branch point singularities on the real axis and this is indicative 
of an imaginary part which vanishes to all orders in the uj'^R/rn^ expansion. In fact. 
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for low frequencies, lmn{uj) is dominated by the closest singularity to the origin, 
leading to the universal behaviour 



271171 

Imnj(co') ~ exp j — - . (4.5) 



The high-frequency limit is again given by (O), where this time the Ci are complex 
but still with Re Cj > 0, so that all polarizations for both Type I and Type II have 
phase velocities that approach c at high frequencies from the superluminal side.^ 



5. Micro- Causality and the Kramers-Kronig Relation 



These results on the analyticity structure of n^"'°°P(ci;) and n{uj) explain why the 
Kramers-Kronig dispersion relation fails to hold for QED in curved spacetime. The 
derivation of relies on the fact that n{uj) is analytic in the upper-half plane. 
As we have seen, however, this is not true in curved spacetime because n{ijj) has 
singularities on the imaginary axis (see fig. (^) and these must be included in the 
contour integral used to derive ( |1 . 1|) . So, for example, in the conformally fiat Type 
I case, the singularities are poles and we have 



duo 



semi-circle ^ 

m{n{oo) -n(0)) +V 



n{uj) - 7rm(0) +P 



■ niuj) 



UJ 



(5.1) 



niuj] 



UJ 



pole contribution 



In this case Imn(a;) = and the principal value integral vanishes, and ( ^.1|) becomes 




Figure 7: The integration contour for f dcu n{uj)/uj used in the proof of the generaHzed KK 
relation for the conformally flat symmetric wave with poles on the imaginary axis. 



^ All this shows clearly that the most important frequency dependence of n{uj) is non-perturbative 
in the parameter uj'^R/m'^. It was therefore not captured by previous effective action approaches [24] , 
which evaluated all orders in a derivative expansion but were restricted to 0{R/ni?). The necessity 
for a non-perturbative technique to determine n(a;) was already noted in refs. [7,24]. 
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Ren(O) -Ren(oo) = ^/ dT {^{1 - J^Res f{t7m/2) , (5.2) 

Jo Jo 

where fix) = (l . ,a ^'^ , a ) ^. Performing the residue sum we find 

•I ^ ' \ smh* X cosh'^ x) x^ o 

aR 

Re n(0) - Re n(oo) = , (5.3) 



which is in perfect agreement with ( [4.1|) and (|4] 

The fact that n{iS) is not analytic in the upper-half plane calls into question the 
validity of micro-causality. Recall that in relativistic QFT, micro-causality, i.e. the 
vanishing of commutators of field operators for spacelike separations, implies that 
the retarded propagator A^et is only non-vanishing in, or, as in the present case of 
a massless quantum, on, the forward null cone. At tree level, this remains true for 
QED in curved spacetime. However, at one loop, the vacuum polarization n]^'°°^ 
contributes to the full propagator: A = A*''*'° — /^tfoe-Qi-ioop^tree -)_... ^ and we must 
check whether n^"'°°P itself only has support in/on the forward cone. 

From our calculation of the on-shell momentum-space vacuum polarization tensor 
n^"^°°P(ci;), we can attempt to determine the dependence of the real-space n^"^°°P on 
the null coordinate v by taking a Fourier transform: 

/oo 
rfcue-^^''ni-^°°P(w) . (5.4) 
-oo 

This is a retarded quantity if the integration contour is taken to avoid singularities 
by veering into the upper-half plane, when f < 0, and the lower half plane, when 
f > 0. For QFT in fiat spacetime, n^"^°°P(co') is analytic in the upper-half plane 
and so when f < one computes the lo integral by completing the contour with a 
semi-circle at infinity in the upper-half plane. Since there are no singularities in the 
upper-half plane, the integral vanishes and consequently Y^~l^°^{y) = for v < 0. 
This is consistent with the fact that the region f < lies outside the forward light 
cone. Hence, in this case micro-causality is preserved as a consequence of analyticity 
in the upper-half plane in frequency space. 

In curved spacetime, on the contrary, H^"^°°p(c(j) is not analytic in the upper-half 
plane and consequently Hjgt°°''(f) may receive contributions from the region f < 
which lies outside the forward light cone. See fig. (§). Indeed, for the conformally fiat 
case where '^~l^°'^{io) has poles on the imaginary axis, we can estimate the behaviour 
of Hj;rP(^) as 



ri-ioop/ ^ /87rm2|t;| 



nr;r^(t^) ~ exp -\\ ^^1^ . (5.5) 
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Figure 8: Including vacuum polarization effects, the photon momentum k may lie outside the 
forward light cone {u > 0, v > 0) of its original null geodesic v = 0. The potential violation of 
micro-causality arises if the retarded propagator is non- vanishing even for v < (the shaded area), 
which lies outside the forward light cone. 

which appears to show a violation of micro-causahty with an exponential dependence 
on a characteristic time/length scale y/R/m^. 

However, in order to be certain that this is a genuine property of the full real- 
space propagator, it is necessary to calculate the vacuum polarization off-shell [25]. 
In our 0{a) on-shell calculation of Iljj, the frequency uj is identified with the light- 
cone component p+, while the component p_ is taken to zero. It remains possible 
that when p_ is small, but non-vanishing, the non-analyticities are shifted into the 
causally safe region Imp+/Imp_ < 0.^ 

It will be especially interesting to explore these issues of superluminal propa- 
gation and causality in spacetimes such as Schwarzschild with both horizons and 
singularities. Although both $oo and vanish for principal null geodesies at an 
event horizon, higher-order terms in the Penrose expansion do play a role and the 
non-vanishing of commutators in the neighbourhood of a horizon could have pro- 
found consequences. The UV-IR effect whereby high frequencies probe the global 
properties of the photon's null geodesic could be especially important in spacetimes 
with singularities. We have implicitly assumed here that the null geodesies in the 
congruence are complete. However, this is no longer true in the presence of singu- 
larities, raising the intriguing possibility that their existence could affect the high 
frequency behaviour of photon propagation in a global rather than purely local way. 
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Dubovsky, Alberto Nicolis, Enrico Trincherini and Giovanni Villadoro for pointing 

Note added: We have now completed a full off-shell calculation of the vacuum polarization 
tensor and find that precisely this behaviour occurs. Full details will be presented elsewhere [26]. 
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sions and Fiorenzo Bastianelli for explaining some details of his work on the world-line 
formalism. This work was supported in part by PPARC grant PP/D507407/1. 

References 

[1] T.J. Hollowood and G.M. Shore, "The refractive index of curved spacetime: the fate 
of causality in QED," arXiv:07072303 [hep-th] 

[2] G. M. Shore, Contemp. Phys. 44 (2003) 503 [arXiv:gr-qc/0304059]. 

[3] G. M. Shore, Nucl. Phys. B 778 (2007) 219 [arXiv:hep-th/0701185]. 

[4] I. T. Drummond and S. J. HathreU, Phys. Rev. D 22 (1980) 343. 

[5] S. W. Hawking and G. F. R. Ellis, "The Large Scale Structure of Spacetime^\ 
Cambridge University Press, 1973. 

[6] S. Liberati, S. Sonego and M. Visser, Annals Phys. 298 (2002) 167 
[arXiv:gr-qc/0107091]. 

[7] G. M. Shore, Causality and Superluminal Light" , in 'Time and Matter', 
Proceedings of the International CoUoquium on the Science of Time, Venice, 
ed. I. Bigi and M. Faessler, World Scientific, Singapore, 2006. [arXiv:gr-qc/0302116]. 

[8] H. A. Kramers, Atti Congr. Intern. Fisici, Como, Nicolo Zanichelh, Bologna, 1927; 
reprinted in H. A. Kramers, Collected Scientific Papers, North-Holland, Amsterdam, 
1956. 

[9] R. de Kronig, Ned. Tyd. Nat. Kunde 9 (1942) 402; Physica 12 (1946) 543. 

[10] S. Weinberg, "The Quantum Theory of Fields^\ Vol /, Cambridge University Press, 
1996. 

[11] R. P. Feynman, Phys. Rev. 80 (1950) 440. 
[12] J. S. Schwinger, Phys. Rev. 82 (1951) 914. 

[13] C. Schubert, Phys. Rept. 355 (2001) 73 [arXiv:hep-th/0101036]. 

[14] F. Bastianelh and A. Zirotti, Nucl. Phys. B 642 (2002) 372 [arXiv:hep-th/0205182]. 

[15] F. Bastianelh, O. Corradini and A. Zirotti, JHEP 0401 (2004) 023 
[arXiv:hep-th/0312064]. 

[16] F. Bastianelh, arXiv:hep-th/0508205. 



- 14 - 



[17] H. Kleinert and A. Chervyakov, Phys. Lett. A 299 (2002) 319 
[arXiv:quant-ph/0206022]. 

[18] H. Kleinert and A. Chervyakov, Phys. Lett. A 308 (2003) 85 
[arXiv:quant-ph/0204067] . 

[19] M. Blau, D. Prank and S. Weiss, Class. Quant. Grav. 23 (2006) 3993 
[arXiv:hep-th/0603109]. 

[20] R. Penrose, "Any space-time has a plane wave as a limit," in: Differential geometry 
and relativity, Reidel and Dordrecht (1976), 271-275. 

[21] M. Blau, "Plane waves and Penrose limits," Lectures given at the 2004 
Saalburg/Wolfersdorf Summer School, 

http : //www . unine . ch/phys/string/Lecturenotes .html 

[22] R. M. Wald, ''General Relativitf , Chicago University Press, 1984. 

[23] G. M. Shore, Nucl. Phys. B 460 (1996) 379 [arXiv:gr-qc/9504041]. 

[24] G. M. Shore, Nucl. Phys. B 633 (2002) 271 [arXiv:gr-qc/0203034]. 

[25] S. Dubovsky, A. Nicolis, E. Trincherini and G. Villadoro, "Microcausality in Curved 
Space-Time," arXiv:0709.1483 [hep-th]. 

[26] T.J. Hollowood and G.M. Shore, in preparation. 



- 15 - 



